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Verfahren zur Erfassung zeitabhangiger Moden 
dynami s che r S y s t erne 

Beschreibung 

Die Erfindung betrifft ein Verfahren zur Erfassung dynamischer 
Systeme, die durch im Zeitverlauf nichtstationare 
Systemparameter charakterisierbar sind, insbesondere ein 
Verfahren zur Segmentierung von Zeitreihen von MeBgroJien 
(Variablen) dynamischer Systeme und zur Identif izierung der 
die Segmente charakterisierenden Systemparameter (Moden) . 

Unter einem dynamischen System wird hier insbesondere jede 
Erscheinung verstanden, deren Zeitverlauf durch eine diskrete 
Abb il dung vom Typ 

x(t+l) = f a(t) (a(t) ) (0.1) 

darstellbar ist. Es werden jedoch auch Systeme mit mehreren 
(z. B. zwei) simultan auf genommenen Zeitreihen x, y gemafi 

y(t+x) = f a(t) (x(t) ) (0.2) 

betrachtet. Dabei bezeichnen a(t) einen Satz 

charakteristischer Systemparameter, x einen Zustand, der im 
allgemeinen einen Vektor in einem mehrdimensionalen 
Zustandsraum bildet, und y einen zeitlich verschobenen 
Zustand. Der Zustandsraum wird von Variablen aufgespannt, die 
z. B. physikalische, chemische, biologische, medizinische, 
geologische, geometrische, numerische und/oder 
prozefitechnische GroBen sein konnen. 

Die Anzahl der Systemvariablen, die zusammen mit der Dynamik f 
das System beschreiben, entspricht der Dimension des 
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Zustandsraumes . Hier werden Systeme betrachtet, deren 

Parameter a ebenfalls zeitlich veranderlich sein konnen . Ein 

gegebenes System mit zeitlich unveranderlichen Parametern a 
wird im folgenden auch als Mode bezeichnet. 

Beobachtbare Oder meflbare Systemvariablen (Mefigroiien) bilden 
erfaflbare Zeitreihen oder Datenstrome, die fur die jeweilige 
Abfolge von Systemmoden charakteristisch sind. Sind die 
Systemparameter innerhalb der Zeitreihen uber bestimmte 
Zeitabschnitte (Segmente) unveranderlich, so kann die 
Zeitreihe entsprechend den jeweils bestehenden Systemmoden 
unterteilt werden (Segmentierung) und jedes Segment einer 
Systemmode zugeordnet werden ( Identif izierung) , 

Viele Vorgange in der Natur wie auch bei technischen 
Anwendungen wiirden sich dann vorhersagen und/oder 
kontrollieren lassen, wenn die ihnen zugrundeliegenden 
dynamischen Prozesse mathematisch modelliert werden konnen. 
Die Analyse und Charakterisierung praktisch gegebener 
dynamischer Systeme wird oft dadurch erschwert, daB sich die 
Systemmoden wahrend der Beobachtung andern. Beispiele hier fur 
sind allmahliche Veranderungen, die sich in Drifts oder Trends 
der Systemparameter widerspiegeln, oder auch spontane oder 
abrupte Veranderungen in der Dynamik komplexer Systeme, etwa 
wenn sich Konf igurationen spontan oder von aufien getrieben 
plotzlich andern. 

Ein Beispiel fur ein betrachtetes System ist die Erzeugung von 
Sprachsignalen im Mund-Rachen-Raum, bei der das System standig 
seine Konf iguration" und - somit seine Mode ^ndert . Es besteht 
ein starkes Interesse daran, die Moden zu erfassen und zu 
- identif izieren, die einem beobachteten Zeitverlauf einer 
Variablen (im Beispiel: Luf tdrucks chwankungen ) zugrundeliegen, 
auch urn verbesserte Vorhersagen uber das betrachtete System 
treffen oder dieses besser kontrollieren zu konnen. 
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Prinzipiell konnen dynamische Systeme anhand vom gemessenen 
Signalen analysiert werden und es sind Reihe von Methoden 
bekannt, aus Zeitreihen Modelle zu gewinnen, die sich fur die 
Vorhersage und eine Kontrolle des Systemverhaltens eignen. So 
ist bekannt, dafi sich der Zustand eines dynamischen Systems 
durch Erfassung der Zeitabhangigkeit beobachtbarer Mefigrofien 
modellieren laiit. Diese Modellierung erfolgt gemafi einem 
ersten Ansatz durch eine Rekonstruktion des Zustandsraumes 
mittels sogenannter Zeit-Verzogerungs-Koordinaten, wie es z. 
B. von N. H. Packard et al . in "Physical Review Letters" (Bd. 
45, 1980, S. 712 ff.) beschrieben ist. Auf der Basis einer 
solchen Rekonstruktion laiit sich dann nur ein einziges 
(globales) Modell f fur die Dynamik finden. Die globale 
Rekonstruktion des Systems ist ferner nachteilig, da bei 
Anwendungen fur vieldimensionale Systeme eine Vielzahl von 
Eingangsgrofien als Randbedingungen vorab bekannt sein mussen 
und/oder aufgrund der hohen Dimensionalitat das System sich 
praktisch nicht mehr schatzen (erfassen, abbilden) lafit 
und/oder ein iibermafiig hoher, nicht praktikabler Rechenauf wand 
entsteht. 

Aufierdem ist dieses Verfahren fur den Fall zeitlich 
veranderlicher Parameter allgemein nicht anwendbar . Die 
Analyse und Modellierung von dynamischen Signalen wird jedoch 
haufig dadurch erschwert, dafi die zugrundeliegenden Systeme 
sich in der Zeit in wesentlichen Parametern verandern. 
Beispiele sind z. B. Signale aus der Medizin, bei denen ein 
Organ wie das Herz oder das Gehirn viele dynamische Moden hat, 
die sich abwechseln, oder etwa den Sprachsignalen, bei denen 

das -er.zeugende- System, .der^Mund-Rachen-Raum^-in-der Zeit — 

of f ensichtlich verschiedene Konf igurationen durchlauft. 

Gemafi einem weiteren Ansatz ist daher aus der Publikation von 
K. Pawelzik, J. Kohlmorgen und K.-R. Muller in "Neural 
Computation" (Bd. 8, 1996, S. 340 ff.) bekannt, Datenstrome 
entsprechend zunachst unbekannten, zeitlich wechselnden 



BEST 




4 



Systemmoden durch Simulierung mit mehreren, untereinander in 
Konkurrenz stehenden Modellen zu segmentieren . Die Modelle 
werden bevorzugt durch neuronale Netzwerke gebildet, die 
jeweils fur eine Dynamik charakteristisch sind und nach 
vorbestimmten Trainingsregeln urn die Beschreibung der 
einzelnen Punkte des Datenstroms konkurrieren . 

Mit dieser Methode gelingt es, eine Zeitreihe in Abschnitte 
quasistationarer Dynamik zu zerlegen und simultan Modelle fur 
diese System-Moden aus den Zeitreihen zu identif izieren . 

Die Segmentierung gemaft K. Pawelzik et al . , zu der unten 
Einzelheiten angegeben werden, erlaubt die Zuordnung von 
Segmenten zu bestimmten Systemdynamiken oder -moden und fuhrt 
zu einer Erfassung des Datenstroms als Vorgang mit einem 
diskreten Schalten (sog. "switching") zwischen den Moden. 
Diese Beschreibung der Parameterdynamik komplexer Systeme 
stellt zwar gegenuber der o. a. globalen Modellierung 
hinsichtlich der Genauigkeit und der Segmentierung 
verschiedener Systemzustande einen Fortschritt dar. Allerdings 
kann der Ubergang zwischen verschiedenen Systemzustanden nicht 
ausreichend beschrieben werden. Es hat sich insbesondere bei 
der Analyse realer Systeme z. B. bei medizinischen Anwendungen 
gezeigt, dali die Segmentierung auf bestimmte Falle mit 
moglichst klaren Modenunterschieden und geringem Rauschen 
beschrankt, allgemein jedoch bei zeitlichen Veranderungen der 
erzeugenden Systeme unzuverlassig ist. 

Solche zeitlichen Veranderungen der erzeugenden Systeme machen 
die beobachtbaren Signals ..nichts.tationar .und. f uhren dazu,. .dafi. . 
sich die Systeme im allgemeinen nicht mehr durch einheitliche 
Modelle beschreiben lassen. Erfolgen diese Veranderungen der 
Systeme plotzlich, so spricht man von Sprungprozessen . 

Die Aufgabe der Erfindung ist es, verbesserte Verfahren zur 
Erfassung der Moden dynamischer Systeme mit nichtstationaren 
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Systemparametern anzugeben, mit denen die Beschrankungen 
herkommlicher Verfahren uberwunden werden konnen und die es 
insbesondere ermoglichen, automatisch mit praktikablem 
Bearbeitungsaufwand und hoher Zuverlassigkeit Segmentierungen 
und Identif izierungen von Zeitreihen mit einer erhohten Zahl 
von Einzelheiten vorzunehmen. 

Diese Aufgabe wird durch das Verfahren mit den Merkmalen gemaB 
den Patentanspruch 1 gelost. Vorteilhafte Ausf iihrungsf ormen 
der Erfindung ergeben sich aus den abhangigen Anspruchen/ 

Die Erfindung basiert auf der Idee, Ubergange zwischen 
verschiedenen Moden eines dynamischen Systems als 
Zwischenmoden des Systems aufzufassen, die paarweise lineare 
Interpolationen der Ausgangs- bzw. Endmoden des Ubergangs 
darstellen. Die betrachteten dynamischen Systeme gehen eher 
allmahlich von einer Mode in die andere liber, statt abrupt 
zwischen Moden zu schalten. Die Erfindung zielt darauf ab, 
solche Ubergange zwischen dynamischen Moden in Signalen und 
die Moden zu identif izieren . 

Es wird daher bei einem Verfahren zur Erfassung der Moden 
dynamischer Systeme z. B. nach einer Schaltsegmentierung einer 
Zeitreihe mindestens einer der Systemvariablen x(t) des 
Systems eine Drif tsegmentierung vorgenommen, bei der in jedem 
Zeitabschnitt, in dem das System von einer ersten Systemmode s 
zu einer zweiten Systemmode Sj iibergeht, eine Folge von 
gemischten Vorhersagemodellen g ± erfafit wird, die durch eine 
lineare, paarweise Uberlagerung der Vorhersagemodelle f i#j der 
.zwei.Systemmoden s^ 3 ...ge.g.eben is.t... . . .. .. 

Gegenstand der Erfindung ist auch eine Vorrichtung zur 
Erfassung eines dynamischen Systems mit einer Vielzahl von 
Moden s t mit jeweils charakteristischen Systemparametern a(t). 
Die Vorrichtung enthalt eine Einrichtung zur Aufnahme einer 
Zeitreihe mindestens einer der Systemvariablen x(t) des 
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Systems/ eine Schaltsegmentierungseinrichtung, die dazu 
eingerichtet ist, in jedem Zeitabschnitt einer vorbestimmten 
Mindestlange fur die Systemvariablen x(t) ein vorbestimmtes 
Vorhersagemodell f ± fur eine entsprechende Systemmode s t zu 
erfassen, und eine Drif tsegmentierungseinrichtung, mit der in 
jedem Zeitabschnitt, in dem das System von einer ersten 
Systemmode s t zu einer zweiten Systemmode ubergeht, eine 
Folge von _ gemischteri Vorhersagemodellen g. erfaflt wird. Die 
erf indungsgemafie Vorrichtung kann ferner Einrichtungen zur 
Einstellung von Interpolations- und Segmentierungsparametern, 
Vergleicherschaltungen zur Verarbeitung von Vorhersagef ehlern 
von Vorhersagemodellen, Anzeige- und Signalisierungsein- 
richtungen und Speichereinrichtungen enthalten. Die 
erf indungsgemafie Vorrichtung kann ein Monitor fur 
physiologische Daten oder physikalische oder chemische 
Prozeliparameter sein . 

Mit der Erfindung steht ein Instrument zur Verfiigung, welches 
ein grofies Anwendungspotential in vielen medizinischen, 
wissenschaf tlichen und technischen Bereichen hat. Mit der 
Segmentierung von Signalen und einer entsprechenden 
Identif izierung der zugrundeliegenden Dynamik eroffnen sich 
neue Moglichkeiten der Vorhersage und der Kontrolle auch 
wesentlich nichtstationarer Systeme. 

Anwendungen der Erfindung haben gezeigt, dafi sich 
kontinuierliche Ubergange zwischen Systemmoden sicher 
identif izieren lieiien und dafi die zugrundeliegenden Dynamiken 
durch die Modelle mit einer Prazision beschreibbar sind, die 

in— vi.e.l.en. Fallen- eine -Vorhersage .des—Systemverhaltens 

ermoglicht. Mit der Erfindung lassen sich in vielen Fallen 
nichtstationarer Prozesse Modelle identif izieren, die sich 
auch fur eine Kontrolle der Prozesse eignen, die ohn'e eine 
Berticksichtigung der Nichtstationaritat nicht moglich waren. 
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Ausfuhrungsformen und weitere Vorteile der -Erfindung werden im 
folgenden unter Bezug auf die beigefiigten Zeichnungen 
beschrieben. Es zeigen: 

Fig. 1: Kurvendarstellungen zur Illustration eines ersten 
Segmentierungsschrittes des erf indungsgemaBen Verfahrens; 

Fig, 2: Kurvendarstellungen zur Illustration eines weiteren 
Segmentierungsschrittes des erf indungsgemaBen Verfahrens 

Fig. 3: Kurvendarstellungen zur Segmentierung von 
Blutregulierungsdaten naeh dem erf indungsgemaBen Verfahren; 
und 

Fig. 4: Kurvendarstellungen zur Segmentierung von EEG-Daten 
nach dem erf indungsgemaBen Verfahren. 

Im folgenden werden zunachst Einzelheiten der Erfindung unter 
Bezug auf die Figuren 1 und 2 und anschlieBend praktische 
Anwendungsbeispiele erlautert. Es ist dem Fachmann 
ersichtlich, daB die Erfindung nicht auf die 

Anwendungsbeispiele beschrankt ist, sondern entsprechend auch 
in anderen Gebieten angewendet werden kann, wie sie 
beispielhaft weiter unten aufgefuhrt sind. 

(1) Erfassung von Drif tiibergangen in nichtstationaren 
Zeitreihen 

Erf indungsgemaB werden nichtstationare Zeitreihen mit einer 

-zweischrittigen Prozedur. ,erf aBt , -bei der . zunachst eine -•- - 

geeignete Modellierung und anschlieBend eine sogenannte 
Drif tsegmentierung erfolgt. Die Modellierung ist dazu 
eingerichtet, in jedem Zeitabschnitt einer vorbestimmten 
Mindestlange fiir jeden Systemparameter ein vorbestimmtes 
Vorhersagemodell fiir eine entsprechende Systemmode zu 
erfassen. Hierzu erfolgt vorzugsweise eine herkommliche 
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Schaltsegmentierung, wie sie beispielsweise aus der 
Publication von K. Pawelzik et al . in "Neural Computation" 
(Bd. 8, 1996, S. 340 ff.) bekannt ist. Die Modellierung kann 
aber auch durch eine andere, in Bezug auf die gewonnene 
Systeminf ormation zur Schaltsegmentierung aquivalente Prozedur 
erfolgen, die an eine konkrete Anwendung z. B. bei bekannten 
reinen Moden oder Randbedingungen angepaiit ist. 

Die Schritte der Schalt- und Drif tsegmentierung werden im 
folgenden naher erlautert. Hinsichtlich der 

Schaltsegmentierung wird hiermit der Inhalt der Publikation 
von K. Pawelzik et al - ausdriicklich vollstandig in die 
vorliegende Beschreibung einbezogen . 

(i) Schritt 1 (Schaltsegmentierung) 

Die Schaltsegmentierung dient der Ermittlung 
charakteristischer Prediktoren, die zur Beschreibung der 
Systemmoden geeignet sind. Die Schaltsegmentierung kann 
entweder an einer Trainings-Zeitreihe oder an der zu 
untersuchenden Zeitreihe durchgefiihrt werden. In beiden Fallen 
konnen die ermittelten Vorhersagemodelle oder Prediktoren fur 
weitere, unbekannte Zeitreihen verwendet werden. 

Es wird ein dynamisches System mit einer endlichen Zahl N 
unterschiedlicher Moden betrachtet. Fur die j-te Mode ist ein 

Wert (allg. : Vektor oder Satz) a 3 (t) eines beobachtbaren 
Systemparameters charakteristisch, der jeweils mit einer 
Funktion f 1(t) (i = 1,...,N) aus einem Satz von N Funktionen f 
modelliert werden -soil— -Es wird die Zeitreihe ~{x t } - =---x j -(-t-)- der 
Systemvariablen betrachtet und zeitabhangig jeweils nach der 
Funktion f i(t) gesucht, fur die {y t } = yj(t) = f i(t} (x.(t)) eine 
neue Zeitreihe vorher zusagender Punkte yj (t) darstellt, die in 
Bezug auf die Systemmoden qualitativ die selben Eigenschaf ten 
hat wie {x t } . Durch den Wechsel der Modellf unktion f im 
Zeitverlauf wird die Schaltsegmentierung gefunden, die die 



BEST AVAILABLE COPY 



9 



Zeitreihe {x t } entsprechend der wechselnden Systemmoden 
unterteilt . 

Die Funktionen f werden als Prediktoren (oder: 
Vorhersagemodelle, Expertenfunktionen) aus einem Satz von 
Netzwerken mit veranderlichen Parametern durch ein geeignetes 
Trainingsprogramm abgeleitet, bei dem simultan sowohl die 
Parameter der Netzwerke als auch die Segmentierung ermittelt 
werden. Der Begriff "Netzwerk" wird hier fur alle moglichen 
geeigneten Modellfunktionen verwendet, also vorzugsweise fur 
neuronale Netzwerke, aber auch z. B. ftlr Polynome oder lineare 
Funktionsapproximationen . Die optimale Wahl eines neuronalen 
Netzwerkes erfolgt in Abhangigkeit von der spezifischen 
Anwendung. Es werden vorzugsweise Netzwerke mit einer 
schnellen Lernf ahigkeit , wie z. B. sogenannte RBF-Netzwerke 
(Radial Basis Function Network) vom Moody-Darken-Typ 
verwendet. 

Das Training erfolgt unter der Voraussetzung, dafl die 
Systemmoden nicht mit jedem Zeitschritt wechseln, sondern eine 
geringere Schaltrate aufweisen, so dafl eine Systemmode liber 
mehrere Zeitschritte erhalten bleibt. Die angenommene Grenze 
der Schaltrate bzw. Anzahl der Zeitschritte, uber die eine 
Systemmode erhalten bleibt, ist zunachst ein freier 
Eingangsparameter und kann anwendungsabhangig in geeigneter 
Weise beispielsweise in Abhangigkeit vorgegebener 
Erf ahrungswerte oder nach einer Parameteranpassungsstrategie 
ausgewahlt werden. Bei der Parameteranpassungsstrategie kann 
vorgesehen sein, fiir die Schaltrate einen Ausgangswert 
- -vorzugeben -und mit diesem . einen Vorhersagef ehler .. (siehe -unten) 
zu ermitteln. Falls die Schaltrate zu hoch oder zu niedrig 
gewahlt ist, fuhrt eine entsprechende Ober- oder 
Unterspezialisierung zu einem zu hohen Vorhersagef ehler . Im 
weiteren Verlauf der Anpassung kann dann die Schaltrate 
optimiert werden, bis der mittlere Vorhersagef ehler unterhalb 
vorbestimmter Grenzen liegt. 
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Das Training erfolgt durch Maximierung der Wahrscheinlichkeit 
W, dafi der Satz von Netzwerken die Zeitreihe {x t } erzeugen 
wiirde. Es handelt sich urn ein Training mit konkurrierendem 
Lernen (engl.: "competitive learning")/ wie es im einzelnen in 
der Abhandlung "Introduction to the theory of neural 
computation" von J. Hertz et al. (Addison-Wesley Publishing 
Company 1991, insb . Kap . 9: "Unsupervised competitive 
learning") beschrieben wird. Aus dieser Abhandlung ist die 
anwendungsabhangige Imp lementie rung eines solchen Trainings 
ableitbar. Die Trainingsregel des konkurrierenden Lernens auf 
der Grundlage des beim Lernen auftretenden Fehlers ist gemafi 
(1) darstellbar. 



Diese Trainingsregel stellt sicher, daB die 
Lerngeschwindigkeit (Verbesserung der Parameter) fur die 
Funktionen f mit geringstem Abstand zum Zielwert y am hochsten 
ist . 

Figur 1 zeigt das Ergebnis der Schaltsegmentierung am Beispiel 
der Analyse einer chaotischen Zeitreihe {x t } mit x t+1 = f (x t ) 
zwischen den vier Moden: 



f x (x) = 4x(l-x) fur x e [0, 1] 
f 2 (x) = f^f^x) ) 

f 3 (x) = 2x fur x g [0, 0.5) bzw. 

-fftxT"- "2 (1-x) tiir-x^e - [0 . 57 1] " 

f 4 (x) = f 3 (f 3 (x) ) 

Fur die ersten 50 Zeitschritte wird mit einem Startwert 
x 0 = 0.5289 zunachst f 1 angewendet. AnschlieBend erfolgt ein 
Obergang (Einzelheiten siehe (ii) ) zur Mode f 2 , die nach 
Schritt 100 bis Schritt 150 stationar wird. Entsprechend wird 




^ j e -P(y-fi) 2 



- (1) 
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ab Schritt 200 bzw. Schritt 300 jeweils fur 50 Schritte die 
Mode f 3 bzw. f 4 eingenommen. Danach erfolgt der Riickubergang 
auf f x . Figur la zeigt einen Ausschnitt (Schritte 300 bis 450) 
des Zeitverlaufs der Zeitreihe {x t } mit x t+1 = f (x t ) . 

Die Seginentierung der ersten 450 Zeitschritte mit 6 
Prediktoren £± f i = 1, . . . , 6 (RBF-Netzwerke vom 
Moody-Darken-Typ) ist in Figur lb dargestellt. Das Training 
ergibt eine Spezialisierung von vier der Prediktoren (6, 2, 4, 
3) jeweils auf eine der o. a. vier Moden. Die stationaren 
Bereiche liegen bei den Intervallen [0, 50] und [400, 450] 
(fj, [100, 150] (f 2 ), [200, 250] (f 3 ) und [300, 350] (f 4 ). Die 
ubrigen zwei Prediktoren (3, 5) haben sich auf die 
Obergangsbereiche zwischen den Moden spezialisiert . Dies zeigt 
den Nachteil der herkommlichen Schaltsegmentierung, bei der im 
Fall von Obergangen der entsprechende Zeitbereich ohne 
adaquate Beschreibung mehrfach unterteilt wird. 

Anstelle des hier beschriebenen sog. "Hard competition"- 
Trainings, bei dem bei einem Trainingsschritt jeweils nur ein 
Vorhersagemodell optimiert wird ("winner takes all") kann auch 
vorgesehen sein, den Grad der Konkurrenz im Rahmen eines sog. 
"Soft competition"-Trainings zu verandern, wie es im einzelnen 
in der Publikation von K. Pawelzik et al . beschrieben ist. 

(ii) Schritt 2 (Drif tsegmentierung) 

Beim zweiten Schritt werden die Ubergange (sog. Drif ten, 
nicht-abruptes ..Ob.exg.ehen, ^gleitendes^Wechseln) zwischen -den _. . . 
Systemmoden beriicksichtigt . Mit der Erfindung wurde als 
wichtige Voraussetzung fur die Drif tsegmentierung 
herausgefunden, dafi der Obergang von einer ersten Systemmode 
direkt in eine zweite Systemmode und nicht uber eine dritte 
Systemmode erfolgt. Das Drif ten zwischen Systemmoden wird 
somit wie folgt durch eine Oberlagerung von (oder paarweise 
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lineare Interpolation zwischen) genau zwei Moden modelliert. 
Dabei treten gemischte, ggf. abgestufte Zwischenmoden auf, die 
jedoch nicht eigene (reine) Systemmoden sind. 

Es wird ein Satz von P reinen Systemmoden, die jeweils durch 
ein Netzwerk k(s), s e P, reprasentiert werden, und ein Satz 
von M gemischten Systemmoden betrachtet, die jeweils durch 
eine lineare Oberlagerung von zwei Netzwerken i(s) und j (s) , 
s eM, reprasentiert werden. Das Modellnetzwerk g 3 fur eine 

gegebene Mode seS, S = P ljM wird gemaB Gleichung (2) 
dargestellt . 

^r^' _ r f k(s)(^t)_ f iir s g P 

g s (x t ) = U(s) f i(s) (x t ) + b(s) f j(s) (x t ) f iir s e M (2) 

In (2) ist it der Vektor (x t/ x t _ T , . . . , x t _ (m _ 1)T ) der 
Zeitverzogerungskoordinaten der Zeitreihe {x t } und f lrj sind 
Prediktoren, die gemaii der o. a. Schaltsegmentierung ermittelt 
wurden. Dabei ist m eine Einbettungsdimension und x der 
Verzogerungsparameter der Einbettung. Die Einbettungsdimension 
ist die Dimension des Phasenraumes, in dem das System 
betrachtet wird und in dem die Modelle operieren. 

Fiir jede gemischte Systemmode sind zwei Parameter a, b 
zusammen mit zwei Net zwerkindizes i, j charakteristisch . Zur 
Vereinf achung des Rechenaufwandes wird die Zahl der gemischten 
Moden begrenzt. Mit 0 < a(s) < 1 und b(s) = 1 - a(s) wird 
eine endliche Zahl von Werten a(s) definiert. Zur weiteren 
Vereinf achung werden "gleiche Abstande zwischen "den "Werten a(s) 
gemali Gleichung (3) gewahlt. 

a r = R I 1 mit r = 1,..., R (3) 

R entspricht der Anzahl der zugelassenen Zwischenmoden und 
wird auch als Auflosung oder Abstufung der Interpolation 
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zwischen den reinen Moden bezeichnet. Die Auflosung R kann 
einen beliebigen Wert annehmen, wird jedoch anwendungsabhangig 
zur Erzielung einer optimalen Systembeschreibung ( insbesondere 
bei stark verrauschten Vorgangen) und praktikabler 
Rechenzeiten insbesondere unter Berticksichtigung der oben 
genannten Schaltrate ausreichend niedrig gewahlt. Bei 
praktischen Anwendungen (siehe unten) kann vorgesehen sein, 
daB die Auflosung R manuell von einem Bediener oder 
automatisch von einem Stellkreis in Abhangigkeit von einem 
vorliegenden Analyseergebnis und einem Vergleich mit einem 
vorgegebenen Schwellwert gewahlt werden . 

Die Gesamtzahl der gemischten Moden betragt bei einer 
gegebenen Auflosung R zwischen jeweils zwei Netzen 
I M | = R * N ' (N-l)/2. Bei dem o. a. Beispiel betragt somit bei 
N = 8 reinen Moden und einer Auflosung R = 32 und Gesamtzahl 
der gemischten Moden |M| = 896. Zur Bestimmung der Gesamtzahl 
der Systemmoden kommen noch die 8 reinen Moden hinzu. 

Die Drif tsegmentierung umfafit nun die Suche nach einer 
Segmentierung mit den reinen und gemischten Systemmoden (a, b, 
R) , die in Bezug auf den Vorhersagef ehler der Moden der 
gesamten Zeitreihe optimiert ist. Die Prediktoren werden so 
ausgewahlt, da£ jedem Element der Zeitreihe einer der Moden 
aus der Gesamtzahl der Systemmoden zugeordnet werden kann. Der 
Vorhersagef ehler ist die Abweichung einer 
Prediktorenvorhersage vom tatsachlichen Element der zu 
untersuchenden Zeitreihe. Fur die zu untersuchende Zeitreihe, 
die nicht mehr unbedingt die Trainings-Zeitreihe ist, mit der 

_bei . _ der _.Schalt segmentierung die __angepafiten.-Netzwerke— oder - - . 

Prediktoren ermittelt wurden, wird fur jeden Zeitschritt mit 
jedem der Prediktoren eine Vorhersage ermittelt, woraus sich 
eine zeitabhangige Matrix der Prediktorenvorhersagen ergibt, 
aus der ein mittlerer Vorhersagef ehler fur willklirlich 
gewahlte Segmentierungen ableitbar ist. Die Segmentierungen 
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mit dem geringsten Vorhersagef ehler ist die gesuchte 
Drif t segment ierung. 

Die Suche nach der Segment ierung mit dem geringsten 
Vorhersagef ehler kann mit jeder geeigneten Such- oder 
Iterationstechnik erfolgen. Vorzugsweise wird eine dynamische 
Programmiertechnik gewahlt, die Equivalent zum 

Viterbi-Algorithmus fur HM-Modelle (sog. Hidden Markov Models) 
ist. Einzelheiten hierzu sind beispielsweise in der 
Publikation "A Tutorial on Hidden Markov Models and Selected 
Applications in Speech Recognition" von L. R. Rabiner in 
"Readings in Speech Recognition" (Hrsg. A. Waibel et al . , San 
Mateo, Morgan Kaufmann, 1990, S. 267-296) beschrieben. Die 
Drif tsegmentierung ist im Rahmen der HM-Modelle die 
wahrscheinlichste Modensequenz, die die zu untersuchende 
Zeitreihe erzeugt haben konnte. Als Nebenbedingung wird dabei 
die Moglichkeit der Moden-Anderungen liber die T-Funktion 
(siehe unten) eingeschrankt . 

Das Ziel der Anpassung ist die Angabe einer optimalen Sequenz 
von Netzwerken bzw. linearen Mischungen aus diesen. Eine 
Sequenz ist optimal, wenn die sogenannte Energie- oder 
Cost-Funktion C* der Vorhersage minimiert ist. Die 
Cost-Funktion C* setzt sich aus der Summe der quadratischen 
Fehler der Vorhersage und den Cost-Funktionen der 
Modeniibergange der Sequenz zusammen. Die Ableitung der 
Cost-Funktion C* jeweils zwischen zwei Zeitpunkten t 0 und t max 
erfolgt induktiv, indem zunachst von einer Start-Cost-Funktion 
gemaB Gleichung (4) ausgegangen wird. 



C 3 (t 0 ) = e 3 (t 0 ) (4) 

wobei 

e a (t) = (x t - g a (l^ t . l ) ) 2 (5) 

der quadratischen Fehler der Vorhersage der reinen oder 
gemischten Moden g ist. 
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Fur den Induktionsschritt von t - 1 auf t wird die 
Cost-Funktion gemafl Gleichung (6) ftir alle s e S berechnet: 

C a (t) = e 3 (t) + min {C s (t-1) + T(s, s)}, t=t 0 +l, • . . t ma x (10) 

seS 

Dabei ist T(s, s) die Cost-Funktion des Ubergangs von einer 
Mode § zu einer Mode s . 

Die optimale (minimale) Cost-Funktion C* ist dann: 

C * = *i? { Cs(t -x)} (ID 

Die Funktion T entspricht bei den HM-Modellen den 
Obergangswahrscheinlichkeiten und kann anwendungsabhangig 
geeignet gewahlt werden. Es kann beispielsweise vorgesehen 
sein, abrupte Schaltubergange und gleitendes Driften zwischen 

zwei Netzen zuzulassen und alle anderen Ubergange durch T = <x> 
auszuschlielien - 

Die Drif tsegmentierung ergibt sich aus der ermittelten 
optimalen Sequenz von Netzwerken bzw. linearen Mischungen aus 
diesen, indem die Moden, die C* ergeben, zurtickver f olgt und in 
Abhangigkeit von der Zeit erfaflt werden. 

Im AnschluB an die Drif tsegmentierung kann sich ein 
zusatzlicher Schritt der Reduzierung der Zahl der zur 
Modellierung verwendeten Netzwerke anschlieflen, der im 
einzelnen unten beschrieben wird. 



AbschlieBend findet eine Identif izierung der segmentierten 
Moden statt, indem jedem Prediktor oder Vorhersagemodell die 
zugehorige Systemmode zugeordnet wird. Diese Identif izierung 
erf olgt anwendungsabhangig . 
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Das Ergebnis der Drif tsegmentierung im Fall der oben unter 
Bezug auf Figur 1 erlauterten chaotischen Zeitreihe {x t } mit 
vier Moden wird im folgenden unter Bezug auf Figur 2 
beschrieben. Die Drif tsegmentierung umfaiJt die Suche nach 
einem Verlauf a(t), der einen speziellen Pfad zwischen den 
reinen Moden liefert, fur den der Vorhersagef ehler der 
gesamten Zeitreihe optimiert ist. 

Jewells nach den ersten 50 Zeitschritte mit der Mode gemaB f 1 
erfolgt fur 50 Schritte ein zeitlinearer Ubergang in die Mode 
gemaii f 2 . Der Ubergang ist ein zeitabhangiges Drif ten gemafi den 
Gleichungen (12) : 



Entsprechende Ubergange finden nach dem 150., 250. und 350. 
Schritt fur jeweils 50 Schritte statt. 

In Figur 2 ist die Besetzung der jeweiligen Moden entsprechend 
den ermittelten Netzwerken in Abhangigkeit von der Zeit 
(Zeitschritte [1200, 2400]) dargestellt. Aus 
Ubersichtlichkeitsgrlinden sind die Ubergangs- oder 
Drif tbereiche entsprechend ihren Zeitgrenzen und Ausgangs- 
bzw. Endmoden in Rahmen dargestellt, in denen der jeweilige 
Driftverlauf zwischen den Moden gepunktet ist. Figur 2a zeigt 
bei einer Auflosung R = 32 (siehe Gleichung (3) ) Ubergange wie 
z. B. bei den Zeitschritten 1350 bis 1400 zwischen den Netzen 
2 und 4. Die Ubergange sind wie gemaii Gleichung (12) erwartet 
linear. Bei einer geringeren Auflosung von R = 3 ergibt sich 
die in Figur 2b dargestellte Segmentierung . Abweichend vom 
linearen Driftverlauf sind die gepunkteten Ubergange 
stuf enf ormig. Dennoch ist die Darstellung auch bei der 
geringeren Auflosung eine adequate Beschreibung des 



f(lt t ) = (1 - a(t))f x (lf t ) + a(t)f 2 ("x t ) 



mit 



(12) 



a(t) = t fc b \\ t a = 50, t b = 100 
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dynamischen Verhaltens des Systems, wie der Vergleich der 
zeitlichen Lage der Moden und des Driftens zeigt. 



(2) Anwendungsbeispiele ftir die Erfassung von Drif tubergangen 

(i) Blut zellenregulation im menschlichen Korper 

Die Blutzellenregulation im menschlichen Korper stellt ein 
hochdimensionales chaotisches System dar, das durch die 
Mackey-Glass-Verzogerungs-Dif f erentialgleichung ( 13 ) ( siehe 
auch in der o. a. Abhandlung von J. Hertz et al . ) beschrieben 
werden kann . 

^P- = -0.1x(t) + °- 2x(t - td 1 ' 0 (13) 

dt x ' l+x(t-t d ) 10 v 1 

Erf indungsgemafl konnen Zeitreihen von physiologischen 
Parametern, die fur die Menge der roten Blutkorper * 
charakteristisch sind, anwendungsabhangig segmentiert werden. 
Die Funktionsf ahigkeit der Segmentierung wird im folgenden 
beispielhaft erlautert. 

Bei Vorgabe von zwei Moden A und B, die sich jeweils durch 
verschiedene Verzogerungsparameter t d = 17 bzw. t d = 23 

unterscheiden, erfolgt mit einer Abtast-Zeitschrittgrofie x = 6 
nach 100 Schritten zunachst ein Ubergang von A nach B. Der 
Ubergang dauert 100 Schritte und ist eine Uberlagerung der 
Gleichung (13) mit den beiden Verzogerungsparameter t d wahrend 
der Integration von Gleichung (13) . Die Uberlagerung wird mit 
einem exponentiellen Drif tparameter a (siehe Gleichung (2)) " 
gemafi Gleichung (14) erzeugt. 

a(t) = exp(f§§) , t = 1,..., 100 (14) 

In der Folge wiederholen sich alle 100 Schritte stationare 
Moden A oder B bzw. die jeweiligen Ubergange. Ftir j eden 
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Riick-Ubergang nach einem Drif tubergang wird ein schaltartiges 
Umspringen angenomiaen. Figur 3a zeigt die entsprechende 
Zeitreihe uber 300 Schritte, Die Drif tsegmentierung mit sechs 
Prediktoren auf der Grundlage von RBF-Netzwerken mit jeweils 
4 0 Basis f unktionen, einem Einbettungsparameter m = 6 und dem 
Verzogerungsparameter x = 1 (siehe Gleichung (2) ) ergibt das 
in Figur 3b gezeigte Bild. Es zeigt sich die erwartete 
Segmentierung der Zeitreihe in stationare Moden und 
Drif tubergange . 

Allerdings zeigt sich, dafl sich jeweils 2 Netzwerke auf eine 
Mode spezialisiert haben (2, 3 => Mode A, 5, 6 :=> Mode B) . In 
einer solchen Situation kann es erf indungsgemaii vorgesehen 
sein, den zusatzlichen Schritt der Reduzierung der Zahl der 
zur Modellierung verwendeten Netzwerke vorzusehen. 

Der Reduzierungschritt umfafit eine sequentielle Reduzierung 
der Netzwerkzahl, jeweils verbunden mit einer Bestimmung des 
mittleren Vorhersagef ehlers . Die Reduzierung (Abzug von 
Redundanznetzwerken) wird beendet, falls eine weitere 
Verringerung der Netzwerkzahl eine signifikante Erhohung des 
Vorhersagef ehlers bedeutet. Figur 3c zeigt das Ergebnis einer 
solchen Reduzierung. Der mittlere quadratische Fehler RMSE 
bleibt bei der Reduzierung um ein, zwei, drei und vier 
Netzwerke jeweils konstant, wohingegen ein starker Anstieg fur 
den Fall auftritt, dali nur noch mit einem Netzwerk modelliert 
wird. Dies bedeutet, dafi das System optimal mit einer Zahl von 
Netzwerken modelliert wird, die gleich der Gesamtzahl von 
betrachteten Netzwerken, vermindert um die Zahl von 
"Reduridanz~netzwerk"en f s t . >~ " ------ 

Die adequaten Modellnetze werden dann erhalten, indem der 
RMSE-Wert fur j ede Net z-Kombination mit reduzierter Netzzahl 
berechnet wird. Die Net z-Kombination mit dem geringsten 
RMSE-Wert umfalit die gesuchten Modellnet zwerke oder 
Prediktoren. Figur 3d zeigt die Drif tsegmentierung nach dem 
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Reduzierungsschritt . Demnach beschreiben die ubrig gebliebenen 
Prediktoren 2 und 5 das System vollstandig. 

( ii ) Schlaf datener fas sung 

Eine weitere Anwendung der Erfindung liegt im Bereich der 
Analyse physiologischer Daten, die fur den Ablauf von Schlaf- 
und Wachmoden von Lebewesen charakteristisch sind. Als 
Grundlage fur sich anschlieflende Prozeduren zur Erfassung von 
Schlaf storungen konnen Zeitreihen z. B. von EEG-Daten 
segmentiert werden . 

Figur 4a zeigt im Vergleich die Ergebnisse einer herkommlichen 
Schaltsegmentierung (oben) , einer Drif tsegmentierung (Mitte) 
und einer "manuellen" Segmentierung (unten) eines 
medizinischen Fachmannes (Schlaf forscher) auf der Grundlage 
von Erf ahrungswerten am Beispiel eines Nachmittagsschlaf es 
eines gesunden Menschen. Die Schalt- und Drif tsegmentierungen 
erfolgen mit acht Netzwerken (netl . . . net8) an 
Einkanal-EEG-Daten x(t) (Figur 4b). In Figur 4a sind wie bei 
Figur 2 aus Obersichtlichkeitsgrunden Rahmen gezogen, die bei 
den Driftmoden verdeutlichen, zwischen welchen Netzen 
interpoliert wird. Die gepunktete Linie im Inneren der Rahmen 
zeigt den jeweils tatsachlichen Verlauf . Die manuelle 
Segmentierung basiert auf der Beobachtung von physiologischen 
Signalen (z. B. EEG, EOG, ECG, Puis, Blutdruck, Atmung, 
Augenbewegungen) . Die Moden Wl, W2 bezeichnen zwei Wachmoden 
mit offenen bzw. geschlossenen Augen und die Moden SI, S2 
^.-jeweils Schlaf zustande- "xu.a . "_aind-."art . " beziehen-.sich- auf — .. 
nicht in Betracht gezogene Zustande bzw. Artefakte. 

Die Schaltsegmentierung zeigt ein vergleichsweise 
undif f erenziertes Bild, das nur grob mit den tibrigen 
Beobachtungen konsistent ist. So tritt beispielsweise ein 
Einschlafvorgang in alien drei Fallen bei t « 7000 auf. Die 
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Drif tsegmentierung ergibt jedoch mehrere Drif tubergange, die 
zusatzliche Einzelheiten des Schlafverhaltens darstellen. Der 
"manuell" beobachtete Schlafbeginn zur Zeit t « 4000 wird 
durch einen exponentiellen Drif tubergang von Netz net7 
(Wachmoden-Prediktor) zum Netz net4 (Schlafmoden-Prediktor ) 

reprasentiert. Das Aufwachen setzt bei t * 9000 durch ein 
leichtes Rtickdriften zum Netz net7 ein, das bis zum Erreichen 

des "manuell" ermittelten Aufwachpunktes t « 9500 gehalten 
wird. In dieser Situation erfolgt eine plotzliche Anderung des 
Wichtungsf aktors, so dafl das Netz net7 eine groBere Wichtung 

erhalt. Nach t « 9800 (Augen of fen) besteht eine Mischung der 
beiden Wachmoden-Prediktoren net7 und net2 . 



(iii) Weitere Anwendungen und Vorteile 

Figur 4a zeigt, daft mit dem erf indungsgemaften Verfahren 
automatisch detaillierte Segmentierungen erzielt werden 
konnen, die bislang nur durch die Beobachtung komplexer 
Merkmalsbilder auf der Grundlage breiter Erfahrungen und 
Intuitionen zuganglich waren. Dieser Vorteil ist nicht nur in 
der Medizin, sondern auch in anderen Gebieten nutzbar, in 
denen grofie Datenmengen bei der Beschreibung komplexer 
dynamischer Systeme anf alien. Derartige Gebiete sind die 
physikalische, chemische und/oder biologische 

Verf ahrenstechnik, die Geologie, Meteorologie, Klimatologie, 
die Spracherf assung u. dgl . . 

Erf indungsgemaBe Verfahren habenjdie folgenden Vo^eile^ Das 
betrachtete System kann hochdimensional sein (10 oder mehr 
Dimensionen) . Die Erfindung erlaubt eine Reduzierung der 
Komplexitat eines solchen Systems durch Betrachtung von 
niedrigerdimensionalen Moden und von wechselnden Ubergangen 
zwischen diesen. Der Einsatz von Vorhersagemodellen fur die 
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Segment ierung ist invariant gegen Anderungen der Amplitude 
erfafiter Signale. 

Die Anwendung der Erfindung fur die Vorhersage oder Kontrolle 
eines Systems erfolgt derart, daft zunachst wie oben 
beschrieben aus der vergangenen Beobachtung und der Kenntnis 
der aktuellen Moden der IST-Zustand des Systems erfafit wird, 
der ggf . eine Mischung entsprechend dem Ergebnis der 
Drif tsegmentierung darstellt. Der IST-Zustand entspricht einem 
dynamischen System f. Die Vorhersage bedeutet, daft das System 
f auf den momentanen Zustand x angewendet wird und sich daraus 
die Vorhersage fur den unmittelbar folgenden Zustand y ergibt . 
Die Kontrolle bedeutet, daft aus dem IST-Zustand die Abweichung 
von einem SOLL-Zustand ermittelt und aus der Abweichung eine 
geeignete Regelstrategie abgeleitet wird. 

Der Vorteil der Vorhersage und Kontrolle besteht darin", daft in 
komplexen Systemen (z. B. bei der Erfassung chemischer 
Reaktionen in einem Reaktor) , die ggf. die Messung lediglich 
weniger Meftgroften erlauben, die fur sich allein aufgrund von 
Mehrdeutigkeiten oder wegen systemimmanenten Verzogerungen 
keine direkten Riickschlusse auf den Systemzustand und ggf. 
bestehende gemischte Zustande erlauben, dennoch detaillierte 
Inf ormationen iiber das System abgeleitet werden konnen. So 
kann bei dem Beispiel mit einer chemischen Reaktion aus der 
erf indungsgemafien Erfassung bespielsweise der makroskopischer 
thermodynamischer Zustandsgroiien eine optimale Regelstrategie 
abgeleitet werden, die die Dosierung bestimmter 
Reaktionspartner umfaBt. 
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Patent anspriiche 

1. Verfahren zur Erfassung der Moden eines dynamischen 
Systems mit einer Vielzahl von Moden s ie die jeweils einen Satz 
<x(t) charakteristischer Systemparameter besitzen, wobei eine 
Zeitreihe mindestens einer Systemvariablen x(t) einer 
Modellierung unterzogen wird, die dazu eingerichtet ist, in 
jedem Zeitabschnitt einer vorbestimmten Mindestlange fiir jede 
Systemvariable x(t) ein vorbestimmtes Vorhersagemodell f ± fiir 
eine entsprechende Systemmode s L zu erfassen, 

dadurch gekennzeichnet, dafi 

nach der Modellierung der Zeitreihe eine Drif tsegmentierung 
erfolgt, bei der in jedem Zeitabschnitt, in dem das System von 
einer ersten Systemmode Si zu einer zweiten Systemmode 
iibergeht, eine Folge von gemischten Vorhersagemodellen g ± 
erfafit wird, die durch eine lineare, paarweise Uberlagerung 
der Vorhersagemodelle f ifj der zwei Systemmoden s i#j gegeben ist. 

2. Verfahren gemafi Anspruch 1, bei dem die Modellierung eine 
Schaltsegmentierung 1st . 

3. Verfahren gemaii Anspruch 2, bei dem die 
Schaltsegmentierung durch Simulierung einer 

Trainings-Zeitreihe des Systems oder der zu untersuchenden 
Zeitreihe mit mehreren, untereinander in Konkurrenz stehenden 
Vorhersagemodellen erf olgt . 

4. Verfahren gemaii Anspruch 3, bei dem die Vorhersagemodelle 
durch neuronale Netzwerke Oder andere Modelle zur Schatzung 
von Funktionen gebildet werden, die jeweils fiir eine Mode s 
charakteristisch sind und nach vorbestimmten Trainingsregeln 
um die Beschreibung der einzelnen Elemente der Zeitreihe 
konkurr ieren . 
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5 Verfahren gemali einem der Anspruche 1 bis 4, bei dem die 
Folge gemischter Systemmoden g, aus den Vorhersagemodelle f . A 
und Interpolationsparametern a, b gemafl 
gi = a(s)fi( S )(x) + b(s)f j(s) (x) ermittelt wird. 

6. Verfahren gemaB Anspruch 5, bei dem die 

Interpolationsparameter gemafi 0 < a(s) < 1 und b(s) = 1 - a(s) 
gewahlt sind. 

7. Verfahren gemali Anspruch 6, bei dem die Werte a(s) auf 
eine bestimmte Auf losungszahl R beschrankt und/oder 
aquidistant sind. 

8. Verfahren gemali einem der vorhergehenden Anspruche, bei 
dem die Folge gemischter Vorhersagemodelle g t erfaftt wird, 
indem fur jeden Zeitschritt mit jedem der moglichen 
Vorhersagemodelle jeweils eine Vorhersage ermittelt wird, 
woraus sich eine zeitabhangige Vorhersagematrix ergibt, aus 
der ein mittlerer Vorhersagef ehler fur willkiirlich gewahlte 
Segmentierungen ableitbar ist, wobei die gesuchte Folge 
gemischter Vorhersagemodelle g L die Segmentierung mit dem 
geringsten Vorhersagef ehler bzw. der maximalen 
Wahrscheinlichkeit ist. 

9. Verfahren gemafi Anspruch 8, bei dem die Suche nach der 
Segmentierung mit dem geringsten Vorhersagef ehler mit einer 
dynamische Programmiertechnik erfolgt, die aquivalent zum 
Viterbi-Algorithmus fur Hidden-Markov-Modelle ist, wobei eine 
optimale Sequenz von Vo r he r s agemode 1 1 en unter Verwendung einer 
minimierten Cost-Funktion C* der Vorhersage ermittelt wird und 
die Segmentierung induktiv aus der Sequenz von 
Vorhersagemodelle abgeleitet wird. 

10. Verfahren gemaii einem der vorhergehenden Anspruche, bei 
dem nach der Drif tsegmentierung ein zusatzlicher Schritt der 
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Reduzierung der Zahl der zur Modellierung verwendeten 
Vorhersagemodelle erfolgt, bei dem die Zahl der 
Vorhersagemodelle sequentiell jeweils verbunden mit einer 
Bestimmung des mittleren Vorhersagef ehlers so weit reduziert 
wird, daB eine weitere Verringerung der Zahl der 
Vorhersagemodelle eine Erhohung des Vorhersagef ehlers 
bedeutet. 



11. Verfahren gemafl einem der vorhergehenden Anspriiche, bei 
dem die Zeitreihe mindestens einer der Systemvariablen x(t) 
eine Zeitreihe von physiologischen Parametern umfaBt, die mit 
der Mackey-Glass-Verzogerungs-Dif f erentialgleichung 

^ = -0.1x(t) + 0 

dt v ' l+x(t-t d ) 10 

beschrieben wird. 

12. Verfahren gemaB einem der Anspriiche 1 bis 11, bei dem die 
Zeitreihe mindestens einer der Systemvariablen x(t) eine 
Zeitreihe von physiologischen Parametern umfaBt, die fur den 
Ablauf von Schlaf- und Wachmoden charakteristisch sind. 

13. Verfahren gemaB Anspruch 12, bei dem die physiologischen 
Parameter EEG-Signale umf assen . 

14. Verfahren gemaB einem der Anspriiche 1 bis 10, bei dem die 
Zeitreihe mindestens einer der Systemvariablen x(t) eine 
Zeitreihe von Sprachsignalen umf aBt . 
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Verfahren zur Erfassung zeitabhangiger Moden 
dynami scher Syst erne 

Zusammenf as sung 

Bei einem Verfahren zur Erfassung der Moden eines dynamischen 
Systems mit einer Vielzahl von Moden, die jeweils einen Satz 
a(t) charakteristischer Systemparameter besitzen, wird eine 
Zeitreihe mindestens einer Systemvariablen x(t) einer 
Modellierung wie z. B. einer Schaltsegmentierung unterzogen, 
die dazu eingerichtet ist, in jedem Zeitabschnitt einer 
vorbestimmten Mindestlange fur jede Systemvariable x(t) ein 
vorbestimmtes Vorhersagemodell wie z. B. ein neuronales 
Netzwerk fur eine entsprechende Systemmode zu erfassen, wobei 
nach der Modellierung der Zeitreihe eine Drif tsegmentierung 

_erfo.lg_t/_bei- der___in _ j edem__Z.eitabschnit.t in-_dem das System, ^von 

einer ersten Systemmode zu einer zweiten Systemmode ubergeht, 
eine Folge von gemischten Vorhersagemodellen erfaBt wird, die 
durch eine lineare, paarweise Oberlagerung der 
Vorhersagemodelle der zwei Systemmoden gegeben ist. 
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